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Abstract. We study global aspects of complete, non-singular asymptotically locally AdS space- 
times solving the vacuum Einstein equations whose conformal infinity is an arbitrary globally sta- 
tionary spacetime. It is proved that any such solution which is asymptotically stationary to the 
past and future is itself globally stationary. 

This gives certain rigidity or uniqueness results for exact AdS and related spacetimes. 

O 
O 

<N 

\ 1. Introduction 

Consider geodesically complete, asymptotically simple solutions of the vacuum Einstein equations 
with negative cosmological constant A < in (n + 1) dimensions. Up to rescaling, these are given by 
complete, (non-singular), metrics g, defined on manifolds of the form M n+1 = R x £, and satisfying 
CO ■ the Einstein equations 

> 

CO ' (1.1) Ricg = -ng. 

The metric g has a conformal completion, at least C 2 , with conformal boundary (X, [7]), where 7 
is a complete Lorentz metric on X. Topologically, conformal infinity X is of the form R x <9£. 
The canonical example is the (exact) anti-de Sitter spacetime gAdSi which may be represented 
<^ . globally in static form as 



X 

5-H 



(1.2) gAdS = ~~ cosh 2 r dt 2 + dr 2 + sinh 2 r g S n- 



1(1): 



<D . . . 

infinity X = R x S n , with boundary metric 70 = —dt + g^n-im the Einstein static cylinder. 



where S l s n - 1 (l) * s the round metric of radius 1 on the sphere S n 1 . Here M = RxR n , with conformal 



Asymptotically simple spacetimes approximate the metric gAdS locally on approach to I, and so 
are often also called asymptotically locally AdS spacetimes. 

It is generally believed that anti-de Sitter spacetime should have an infinite dimensional space 
of dynamical perturbations, i.e. time-dependent, complete vacuum solutions (1.1), which have the 
same conformal infinity (X, 70) as exact AdS and which are globally close to gAdS- I n other words, 
gAdS is dynamically stable. This is certainly the case at the linearized level; one may globally solve 
the linearized Einstein equations at gAdSi with zero boundary data on X and arbitrary smooth 
Cauchy data on S. These linearized solutions, (or normalizable modes), remain uniformly bounded 
in time, cf. [1] for a detailed treatment. 

Such dynamical or global stability results are well-known in the context of spacetimes with zero 
cosmological constant, A = 0. Thus, the celebrated work of Christodoulou-Klainerman [2] shows 
that there exist global non-singular perturbations of Minkowski spacetime, (in 3+1 dimensions), 
which tend to the flat Minkowski spacetime as t — > ±00; see also [3] and the more recent work [4], 
as well as [5], which gives the existence of non-singular asymptotically simple global perturbations. 
There are also analogues of such stability results in the context of cosmological spacetimes, (in the 
expanding direction), where E is compact without boundary, cf. [6], [7]. 
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Similarly, in the context where A > 0, Friedrich [8] has proved that exact de Sitter spacetime is 
globally stable in a natural sense in 3 + 1 dimensions; the same result holds in fact at least in all 
even dimensions, cf. [9]. 

In this paper, we discuss some aspects of the global dynamics of complete, asymptotically simple 
solutions of the Einstein equations when A < 0. To describe the main results, let 7 be a fixed 
but arbitrary globally stationary metric on X ~ R x <9£, with <9X compact. Let £ = £(7) be the 
space of all geodesically complete, asymptotically simple solutions of the Einstein equations (1.1) 
with conformal infinity (X, 7); (this definition will be made more precise in §2). Choose a fixed 
diffeomorphism M ~ R x S, (i.e. spacetime decomposition), and let C denote the corresponding 
space of solutions of the constraint equations on (T,,g,K) given by solutions in £. In particular, a 
metric g G E gives rise to a curve (gt, K t ), t G (—00, 00) in C. 

Theorem 1.1. Suppose g G £ is weakly asymptotically stationary in the sense that there exist times 
U — > 00 and ti — > —00 such that (g^^K^) converges to a globally stationary solution of the vac- 
uum constraint equations, as i — > 00. Suppose further that, modulo infinitesimal diffeomorphisms, 
solutions of the linearized Einstein equations at g satisfy the unique continuation property at X. 
Then (M n+1 ,g) is globally stationary. 

As explained in detail in §2, the unique continuation property at X means that solutions of the 
linearized Einstein equations are uniquely determined, up to infinitesimal diffeomorphisms, by their 
Cauchy data, {not just the boundary data), at X. Equivalently, they are uniquely determined, up 
to infinitesimal diffeomorphism, by their formal series expansion at X. This is automatically the 
case for solutions of the linearized equations which are analytic, (in the polyhomogeneous sense), 
at X. The unique continuation property was proved to hold in general for Riemannian-Einstein 
metrics in [10], and we conjecture that it also always holds for Lorentzian-Einstein metrics. There 
certainly appears to be no physical reasons to doubt the validity of this property. 

The time evolution of metrics g G E gives rise to a flow, (i.e. dynamical system), on the constraint 
space C. Stationary solutions in £ then correspond to fixed points of this flow, (for a suitable choice 
of spacetime foliation of M). Roughly speaking, Theorem 1.1 then states that there are no orbits 
of the flow on C which are weakly asymptotic to a fixed point at both ends, except fixed orbits. 

It is possible that (X, 7) has more than one timelike Killing field, (modulo constants). This will 
occur when there are non-trivial spatial Killing fields on (X, 7). Theorem 1.1 applies to each Killing 
field and leads easily to the following corollary. 

Corollary 1.2. Under the assumptions in Theorem 1.1, suppose (X, 7) is the Einstein static cylin- 
der and that (M n+1 ,g) is weakly asymptotic to the exact AdS spacetime to the infinite past and 
infinite future. Then (M n+1 ,g) is globally isometric to the exact AdS spacetime. 

Similar rigidity or uniqueness results hold for other spacetimes which, for example, have sufficient 
symmetry at conformal infinity (X, 7). Thus, the Horowitz- Myers AdS soliton [11] is unique in 
the sense above, cf. §3. The same applies to the AdS soliton metric of Copsey-Horowitz [12]. 
Previously, the uniqueness of the Horowitz-Myers soliton metric was only known within the much 
more restrictive class of static solutions of the Einstein equations, cf. [13], [14]. 

A more general result on the extension of isometries from the boundary to the bulk is proved in 
[10] for Riemannian (or Euclidean) Einstein metrics (1.1), by completely different methods. 

The proofs of the results above are conceptually very simple. They follow from basic conservation 
properties of the holographic stress-energy tensor and holographic mass arising in the AdS/CFT 
correspondence, [15], [16], [17], [18], [19], together with a basic identity discussed in §3, (cf. (3.1)). 
For simplicity, we have restricted the analysis to vacuum solutions of the Einstein equations (1.1) 
with negative cosmological constant. Of course the results above also hold in the presence of 
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matter terms which arise from a vacuum solution in higher dimensions via Kaluza-Klein reduction. 
However, in general the presence of matter terms changes the conservation properties of the stress- 
energy tensor [17]. In any case, we hope to discuss the situation with matter terms elsewhere. 

As stated, these results do not hold for spacetimes containing black holes, or for any spacetimes 
containing singularities which propagate to X. As a specific example, it is well-known that the 
AdS-Kerr spacetime has conformal infinity X given by a finite-time region / x S n - X of the Einstein 
cylinder R x S n ~ x . Not all symmetries of S n ~ l extend to symmetries of the AdS-Kerr solution, 
so there are timelike Killing fields on X which do not extend to Killing fields on the Kerr-AdS 
solution. More generally, smooth global initial data on a space-like slice £ may evolve to the 
future (or past) and eventually form a black hole, for example the AdS Kerr metric. Such solutions 
will be asymptotically stationary to the future (or past) near X within their domain of outer 
communication, but are not themselves globally stationary; see also Remark 3.1. 

The contents of the paper are briefly as follows. In §2, we present some basic background material 
and elementary results needed to establish the main results. These results are then proved in §3. 
In §4, we conclude with further discussion and interpretation of the results. 

I would like to thank Piotr Chrusciel, Gary Gibbons, Gregory Jones for discussions and corre- 
spondence on issues related to this work. Particular thanks to Kostas Skenderis and Bob Wald for 
numerous comments and criticisms on an earlier version of this work. This project began during 
a stay at the Newton Institute, Cambridge in Fall 05, and I am grateful to the Institute for its 
hospitality. 

2. Background Material. 

Let X be a compact n-manifold with boundary <9X and M = KxE, where M. represents the time 
direction. We consider geodesically complete solutions (M, g) of the vacuum Einstein equations 
(1.1) which are asymptotically simple. Asymptotic simplicity is equivalent to the existence of a 
(reasonably smooth) conformal completion. Thus, let p be a defining function for the boundary 
dM = I = Ix 9S, i.e. p is a coordinate function for dM in M, vanishing on dM. One then 
requires that the conformally equivalent (unphysical) metric 

(2.1) ~g = fg 

extends at least C 2 to X; (a stronger smoothness condition will be required below). It is easy to see 
that this implies that the sectional curvatures K of g satisfy \K + 1| = O^p 2 ), so that the metric 
g locally approaches the AdS metric to order 0(f) 2 ) near X; for this reason, asymptotically simple 
spacetimes are also often called asymptotically locally AdS. The space X is called conformal infinity. 
Let 

(2-2) j = g\j 

be the metric induced on X. It is well-known that causality arguments imply that 7 is a Lorentz 
metric on X. Different choices of defining function p lead to conformally equivalent metrics 7 on X, 
so that only the conformal structure [7] on X is uniquely determined by (M, g) . 

If (M,g) is asymptotically simple, then it is standard, (and easily seen), that each choice of 
boundary metric 7 £ [7] determines a unique geodesic defining function p, for which the integral 
curves of Vp in the compactified metric 

(2.3) g = p 2 g 

are (spacelike) geodesies orthogonal to X. In the following, we work only with such geodesic 
compactifications. (The integral curves of Vlogp are also geodesies with respect to the metric g). 
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The metric 5 splits in the p-direction, so that 

(2.4) g = p- 2 (dp 2 + g p ), 

where g p is a curve of Lorentz metrics on the level sets S(p) of p. As p — > 0, g p — > 7 = 5( ). The 
Feffer man- Graham expansion [20] is the expansion of the curve g p in a Taylor-type series in p. The 
exact form of the expansion depends on whether n is even or odd. If n is odd, then one has 

(2-5) g p ~ 5(o) + p 2 5(2) + • • • + P n_1 5(n-i) + P n 5(n) + ■ ■ ■ 

while if n is even, 

(2-6) 5p ~ 9(0) + P 2 5(2) + • • • + P n ~ l 9{n-l) + P n log P H + P n 5(n) + " " " 

Below order n, the expansions are even in powers of p, and all coefficients 5(2fc), 2/c < n, as well as 
the coefficient H in (2.6) are determined by the boundary metric 7 = 5( ) and its derivatives up to 
order 2k, (respectively n); thus they do not depend on the particular bulk Einstein metric g. The 
series (2.5) is a formal power series, in powers of p, while the series (2.6) is a series in powers of p 
and log p, (i.e. a poly homogeneous series). 

The coefficient 5( n ) is formally undetermined; thus it is not determined by 7, and depends on the 
bulk metric g. All the remaining (higher order) coefficients in the series (2.5)-(2.6) are completely 
determined by the data 

(2-7) (s( ), </(„)), 

so that these terms determine the formal expansion of the metric 5, and hence g, near X. Note 
that the terms in the expansion (2.5)-(2.6) depend on the choice of boundary metric 7 € [7]. A 
conformal change of 7 will cause a change in the geodesic defining function p, and hence a change 
in the coefficients. Transformation formulae for these coefficients are given in [17], [21] for instance. 

The boundary metric gr \ is (formally) free, in the sense that the Einstein equations in a neigh- 
borhood of I impose no conditions or constraints on g^ . Similarly, the transverse-traceless part of 
<7( n ) is formally free; however, the Einstein equations do impose constraints on the divergence and 
trace of 5( n ) . Thus 

(2.8) <%«)=?>), and trg {n) = s (n) , 

where the divergence 5 and trace are taken with respect to gr \. The terms ri n ) and sua may be 
explicitly computed from the boundary metric 5(0) and its derivatives up to order n, cf. [17]. When 
n is odd, one has r^ n ) = S( n ) = 0. We will call the equations (2.8) the constraint equations on X. 
They arise from the Gauss and Gauss-Codazzi equations on 1 in the geodesic gauge (2.4). 

Throughout the following, we assume that the metric g is asymptotically simple in the sense 
that the expansion (2.5) or (2.6) exists to order n, so that 



(2-9) g p = 5(0) + P 2 9(2) + ■■■ + p" _1 5(n-i) + P n 9(n) + o{p n ), 

or, 

(2.10) g p = 5(0) + P 2 5( 2 ) + • • • + P n_1 5(n-i) + p n log pH + p n g (n) + o(p"), 

where o(p n )/p n — > as p — > 0. 



We point out that if the free data (5(o)>5(n)) are real-analytic and satisfy the constraints (2.8), 
then a result of Kichenassamy [22] shows that the (formal) series expansion (2.5) or (2.6) converges 
to an actual solution 5 of the Einstein equations near X. Of course, since the series is uniquely 
determined by the data (2.7), this solution is unique. (It is mathematically an open question 
whether there could exist other solutions with the same data (2.7) which don't have convergent 
expansions; this seems very unlikely however). 
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Throughout most of the paper, we only consider boundary metrics 7 = <7( ) which are globally 
stationary on T. Thus, there is a complete, timelike Killing field Z, generating a free isometric 
]R-action on (1, 7), so that 

(2.11) £zl = 0. 

Let 7r : I — > S be the projection onto the orbit space of the R-action. The M-bundle ir is trivial, 
M ~ R x S, and with respect to a fixed trivialization determined by a global time function i on 1, 
the metric 7 may be written globally in the form 

(2.12) -y = -N 2 (dt + e) 2 + ir*g s , 

where 9 is a connection 1-form on the bundle n, N 2 = —j(Z, Z) > and gs in the Riemannian 
metric induced on the orbit space S by 7. 

We require that 7 is C n+1,a smooth, but otherwise impose no other conditions on 7; it need 
not satisfy any equations or have any other symmetries. If 7 does have other symmetries, i.e. the 
metric gs also admits Killing fields, the timelike Killing field Z will not be unique. The results to 
follow hold for any fixed choice of Z. 

It follows that if g is an asymptotically simple solution of the Einstein equations with boundary 
metric 7, then the determined coefficients in (2.5)-(2.6), i.e. gi^ with k < n — 1, and the logarithmic 
coefficient TC, are also invariant under the isometric action generated by Z, 

(2.13) C z g (k) =0, CzH = 0. 

Definition 2.1. The space 8 = £ 7 is the space of all geodesically complete, asymptotically simple 
solutions of the vacuum Einstein equations which have the fixed stationary metric 7 in (2.12) as 
boundary metric and for which the expansions (2.9)-(2.10) hold to order n+ 1 + a, a G (0, 1). One 
may define a natural C n+1 ' a poly homogeneous topology on S . 

Let t be a fixed smooth global time function on M, which restricts to the time function t above 
on X. This gives a trivialization M = R x S, with fibers S^ given by the level sets of t. With 
respect to this foliation S^, the metric g may be written in local coordinates (t,x l ) as 

(2.14) g = -u 2 dt 2 + g ij (dx i + Cdt)(dx j + tfdt), 

where (TV, £) is the lapse-shift vector with respect to dt- Since the boundary metric of g with respect 
to p is 7, one has up — > TV, for N as in (2.12), as p — > 0. 

Let C be the space of solutions of the constraint equations on a given S, say So as above, induced 
by the global metrics g € £. Thus, an element in C is a triple (T,, go, K), where go is a complete, 
conformally compact metric on S, and K is a symmetric bilinear form, (the 2nd fundamental form), 
with K = go + 0{p 2 ). The data (go,K) satisfy the vacuum constraint equations 

(2.15) R- \K\ 2 + k 2 =n(n-l), 

SK + dK = 0, 

where n = trK is the mean curvature. 

Note that C does not consist of all solutions of the constraint equations (2.15); many solutions 
of (2.15) will only give rise to local-in-time solutions of the vacuum equations (1.1). The space C 
may be given a C n+1,a x C n,a polyhomogeneous topology induced from the topology of 6. 

Given the spacelike foliation S t above, let gt be the metric induced by g G £ on S t . Choosing 
the zero-shift gauge gives a diffeomorphism <pt '■ S = So — > S^, and we will let gt also denote the 
induced metric <ftt(9t) on S. Thus, the Einstein flow from S to any S< defines a flow, denoted 

(2.16) g - ft, K -> X t 
on the constraint space C. 
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Next we discuss briefly the linearized Einstein equations. The linearization of the vacuum equa- 
tions (1.1) at g € £ is given by 

(2.17) D*Dh-2Rh + 5*5h = 0. 

Here h is a symmetric bilinear form, h = h — \trhg, D*D = trD 2 is the wave operator and R is 
the curvature tensor acting on symmetric bilinear forms; all metric quantities in (2.17) are with 
respect to the background metric g. Altering h by infinitesimal diffeomorphisms h — > h + S*X, it is 
well-known that one can solve (2.17) in the transverse gauge 5h = 0, by solving the coupled system 

(2.18) D*D(h + S*X) - 2R(h + 5* X) = 0, 

(2.19) 5{h + 5*X) = 0, 

in the variables (h,X), cf. [23]. 

Consider the initial boundary value problem for (2.18): 

(2.20) D*Dh - 2Rh = 0, % = ho, V dt h = hi, and h\j = h {0) , 

where the boundary data fy ) e C°°(Z) and the initial data ho, hi are C°° poly homogeneous on 
X = So up to the boundary. We assume ho, hi and hn^ match in a smooth poly homogeneous 
sense at the corner <9£. In this generality, it is not known if there is a global C°° poly homogeneous 
solution h of (2.20) defined on all of M. 1 However, it suffices for our purposes, (namely regarding 
the unique continuation property), that one has a smooth solution in the interior of the domain 
of dependence D(T,) C M, which is C n+1 ' a polyhomogeneous smooth at the boundary <9£ C I. 
The existence of such solutions follows from standard energy estimates, cf. [24] for instance. Such 
energy estimates are carried out in the Sobolev spaces H s ; note that Sobolev embedding gives 
(jn+i,a (- jjs^ £ Qr g y | n _|_2. in fact, it suffices for our purposes to know that there exist solutions 
to the linearization of the constraint equations (2.15) which are C n+1,a polyhomogeneous at <9£; 
this follows from the work of [25]. 

A standard identity (Weitzenbock formula) gives 55* X = \D*DX — RicX, so that (2.19) is 
equivalent to the vector system 

(2.21) \D*DX - RicX = -5h, 

which is of the same form as (2.20). Thus, as above, one may solve the initial boundary value 
problem for (2.21) within -D(S). 

The symmetric bilinear forms h satisfying 

(2.22) D*Dh - 2Rh = 0, 

5h = 0, 

which satisfy h/ ) = on J may be viewed as defining the formal tangent space T g £ to £, (modulo 
diffeomorphisms). However, this is formal; no claim is made that £ is a manifold, corresponding 
to the linearization stability of the Einstein equations within the space £. Moreover, solutions of 
(2.22) with 

h(o) + on J, 

are certainly not in the formal tangent space T£, since these deformations don't preserve the 
boundary metric. 

We now define precisely the unique continuation property (iii) of Theorem 1.1. This is the 
statement that, for g G £ , any solution h of the linearized Einstein equations (2.22) which vanishes 

^The logarithmic terms in the expansion of h at I may propagate into the bulk of M and could, apriori, lead to 
singularities of the solution in the bulk. I am grateful to P. Chrusciel for discussions on this point. 
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to infinite order at X is necessarily zero. For solutions h of (2.22) which are C°° poly homogeneous 
up to X, this is clearly equivalent to the statement that if h has zero Cauchy data on X, i.e. if 

(2.23) h( ) = and h (n) = 0, 
then 

(2.24) h = 0. 

This follows from the properties of expansions the (2.5)-(2.6); the conditions (2.23) imply that 
the formal series solution of (2.22) vanishes. Note that it suffices to have the unique continuation 
property at a cut dH of X, i.e. within the domain of dependence £>(£). For if h vanishes to 1st 
order on £ and = 0, then h = 0, by uniqueness of solutions to the initial boundary value 
problem (2.22). 

Next we return to the expansions (2.5)-(2.6). The undetermined term is closely related to 
the stress-energy tensor on X, and is an important feature of the AdS/CFT correspondence. Thus, 
as shown in [16], [18], there is a symmetric bilinear form r( n ), determined by the boundary metric 
7 and its derivatives up to order n, such that the form 

(2-25) T (n) = # (n) - r (n) , 

is divergence-free with determined trace, i.e. 

(2.26) 5 7 T (n) = 0, tr 7 r (n ) = a (n) . 

The term T( Tt ) is obtained by a covariant (intrinsic) renormalization of the Brown- York quasi-local 
stress-energy tensor. Via the AdS/CFT correspondence, T( n ) corresponds to the expectation value 
of the stress-energy tensor of the QFT dual to (M, g) . The term a( n ) is proportional to the conformal 
anomaly [26], and is determined by the boundary metric 7. When n is odd, rr n \ = 0, so TY n ) = gt n \ 
and a( ra ) = 0. 

It is important to note that the construction of T( n ) is background- independent; there is no 
normalization with respect to a background "standard" solution. Indeed, there are no such standard 
background solutions to which (M, g) can be compared in the generality of the current discussion. 
For instance, apriori, there may not be any stationary metric on M with conformal infinity (X, 7). 
A more recent and efficient construction of the stress-energy tensor 77 n ) is given in [18], [19], cf. also 
[27], [28]. 

Given g € £, the (holographic) mass, cf. [16], [17], of the cut dT, t C X is defined by 

(2.27) m dSt = T( n) (Z,v)dVy, 

JdT, t 

where v is the future unit normal and Z is the (future-oriented) Killing field on X. By its definition, 
this mass is independent of spacelike hypersurfaces S t C (M, g) giving the same cut dT> t at 2. Both 
the holographic mass m and the stress-energy tensor T( n ) depend on the choice of boundary metric, 
(or equivalently on the choice of defining function p). However, as noted above, the boundary 
metric 7 is chosen to be the fixed stationary metric (2.12). 

Since 7 is stationary, a standard application of (2.26) and the divergence theorem implies that 

(2.28) m = m dSt , 

is independent of the cut <9£ t , (and hence t), so that the mass depends only on the solution 
(M,g) £ £, (given the fixed choice of boundary metric), and the choice of timelike Killing field Z 
on (X, 7). In other words, the mass is conserved. If 7 has a larger space of timelike Killing fields, one 
may consider the mass (2.27) defined with respect to each choice, or take the "canonical" Killing 
field with zero angular velocity, cf. [19]. For the rest of the paper, Z denotes any fixed timelike 
Killing field. 
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There have been numerous definitions of mass and other conserved quantities for asymptotically 
AdS and asymptotically locally AdS spacetimes, cf. [29], [30], [31], for example; cf. [32] for an 
overview. In the generality considered here, the holographic mass in the only suitable definition, 
again since it is background- independent and covariant or intrinsic to the metric (M,g), given a 
fixed choice of boundary metric. Many of the various definitions of mass have recently been shown 
to be equivalent to the holographic mass, cf. [33] and in particular [19] for a very clear analysis. 

Next, consider variations of the mass with respect to variations of the metric g £ £. Suppose 
then that h is an infinitesimal Einstein deformation, with induced boundary variation fy )- It is 
convenient to put the variation h in the geodesic gauge for a fixed geodesic defining function p, so 
that h respects the splitting (2.4), i.e. ho a = 0. This can always be accomplished by a suitable 
gauge transformation, (i.e. infinitesimal diffeomorphism) . Since the curve g u = g + uh is Einstein 
to 1st order in u, the constraint equations (2.26) give 

(2.29) 5'(T {n) ) + 5a (n) = 0, 

where 5' is the variation of the divergence 5 in the direction of the variation of the boundary metric 
7 and <T( n ) = = /i( n ) — r'^ n y Suppose the boundary metric 7 is kept fixed, i.e. /i( ) = 0, so that, 
formally, h G T g 8. Then 

(2.30) 5a (n) = 0. 
By definition, one has 

(2.31) dm d x t (h) = ^^\ u=0 = [ -^-(T (n) ) g+uh (Z, u)dV = [ (j {n) (Z,v)dV, 

again since the boundary metric is fixed. Thus, via (2.30) and the divergence theorem, drriQ^ is 
independent of the cut <9E, i.e. is conserved, for any fixed h vanishing on 1. So for such variations, 
both the mass as well as its variation, are conserved. 

The following discussion will be important in the proof of Theorem 1.1 in §3. For an arbitrary 
complete Lorenztian boundary metric 7 on I = M x <9X, let T be the space of (smooth) symmetric 
bilinear forms on 2 which satisfy the constraint equations (2.26) with respect to 7; when n is odd, 
these are the transverse-traceless (TT) forms. Thus T is naturally an affine bundle 

7T : T -> Met(l) 

over the space of Lorentz metrics on J. An element (7, r), r <G T 7 , in T then defines a unique formal 
series solution of the Einstein equations (1.1), defined near I. This follows immediately from the 
discussion following (2.5)-(2.6). If the pair (7, r) are real-analytic, the formal series converges to 
an actual solution g of the Einstein equations (1.1), again defined in a neighborhood of J. 

Proposition 2.2. At any (7, r) G T, the map it is a submersion, i.e. its derivative is surjective, 
and so ir is locally surjective. 

Proof: Given (7, r) £ T, one needs to show that for any variation 7' = /i( ) of 7, there exist 
solutions t' of the linearized constraint equations (2.26), i.e. 

(2.32) S't + St' = 0, 

tr' t + trr' = a'^ . 

Let a = t', so that it suffices to solve 



(2.33) 



5a = 4>i, tra = <p 2 , 
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for arbitary (pi,(p2- Consider for instance a of the form a = 5*V + /7, where V is a vector field. 
Then (2.33) becomes the system 55*V - df = <f> u tr(5*V) + nf = <p 2 , so that 

5S*V + ±dtr6*V = </>!- ±d<f> 2 . 

One has 55* V = \D* DV ' — RicV ', while tr 5*V = (^ — 1)5V. Hence, the equation above is equivalent 
to 

(2.34) \D*DV + {\- l)d5V - RicV = fa- ±# 2 . 

This is a linear hyperbolic system for V on I with dT, a compact Cauchy surface. It is standard that 
the Cauchy problem for (2.34) has global solutions V on (1, 7) with arbitrary initial data. Given 
V, one may then solve the trace equation above to obtain / and the resulting pair a satisfying 
(2.32). 

■ 

Just as before with pairs (7, r) satisfying the constraint equations (2.26), the space of pairs 
(h(fy,a), for a satisfying (2.32), corresponds exactly to the space of formal series solutions of 
the linearized Einstein equations defined near I, (in the geodesic gauge). Again if (h^,a) are 
real-analytic, the formal series converges to an actual solution of the linearized Einstein equations 
defined in a neighborhood of 2. 

For a given fy )> the space of solutions a of (2.32) is an affine space ^ht y Observe that the 
space of global solutions, (or solutions defined in £>(E)), of the Einstein equations linearized at g 
with induced boundary variation fy ) i s a l so an affine space £/ft (0) , parameterized by the Cauchy 
data /to, h\ on a spacelike hypersurface E C (M,g). Clearly, 

We do not know if the two spaces Gh (0) , ^Vi( ) actually coincide, although there seems to be no 
compelling reason for this to be the case. 

Lemma 2.3. Suppose 7 is globally stationary. Then for any smooth variation htQ\ of compact 
support on 1, there exist solutions a to (2.32) on X which are uniformly bounded, in that 

(2.35) \a\ c i < K, 
where K depends only on (7, fy )) an d is independent oft. 

Proof: We will constuct a specific bounded solution, although it is clear that there will be 
many other possibilities. Let e a , < a < n, be a local orthonormal framing of (X, 7), with 
eo = T = Z/\Z\, so that {e^} are tangent to the orbit space S oil. At a given point p £ S, assume 
that V ei ej{p) = 0, where here, and only here, the covariant derivative is on (S, 75-). Then, at p, 
-5a = (v T a)(T) + (V ej <j)( ei ) = V T <r(T) - ^(VlogiV) + V e ^( ei ), since (V^T) = 0. (The last 
equation follows since T is colinear to a Killing field orthogonal to ej). Let o~ H be the horizontal 
part of a, a 11 = a - a(T) ■ T. 

For simplicity, suppose a(T, e«) = 0, for all i. Then a simple calculation shows that the first 
equation in (2.33) becomes 

(2.36) -(Vt<t(T))(T) = <I>(T), 

(2.37) 5 H (Na H ) = N(j) H , 

where (j) H = <fi — <fi(T)T is the horizontal projection of (j) and 5 H is the divergence operator on the 
orbit space (S,gs). Since /i( ) has compact support, so does <fi. 

9 



The equations (2.36)-(2.37) are uncoupled. The first equation may be solved directly by integra- 
tion along the integral curves of T. Since <j)(T) has compact support, a(T,T) remains uniformly 
bounded in the sense of (2.35). 

For the second equation, let E> 2 (S) be the space of symmetric bilinear forms on S. Then S 2 (S) = 
Im{5 H )®Ker{(6 H )*) : and the second summand corresponds to the space of Killing fields on (S, 75). 
Suppose first the orbit space (S,gs) has no Killing fields, so that the operator 6 : § 2 (S) — > $7 X (5) 
is surjective. Then (2.37) admits (many) solutions, with a H satisfying (2.35). 

If (S, 75) admits Killing fields, then by linearity it suffices to solve 

(2.38) 6a = X t , 

where, for each fixed t, X t is a Killing field on S and X t has compact support on X. Consider first 
the equation 

6a = f(t)X, 

where X is a fixed Killing field on S. Setting a = v(t)T ■ X, a simple computation gives 

-6{v{t)T ■ X) = v'{t)X + V X T + V T X. 

Straightforward computation using the Killing properties shows that V jT + VtX is orthogonal to 
T and the space X (S) of Killing fields on S, and hence V X T + V T A £ Im{6 H ). Setting v' = -f 
and using linearity then shows that the first equation in (2.33) is solvable with bounded a. A 
general Xt as in (2.38) has the form Xt = Yl fi(t)Xi, where X{ is a basis of x(S) and so the result 
for the first equation in (2.33) again follows by linearity. 

To solve the second equation in (2.33), by linearity it suffices to solve 6a = and tra = 4>2- It 
is clear that this system has many C 1 bounded solutions. 

■ 

It is not clear in this generality that there exists a solution h of the linearized Einstein equations, 
with h\x = /i(o) of compact support for which <7( n ) satisfies (2.35). 



3. Proofs of the Results. 

The main tool used in the proof of Theorems 1.1 and 1.2 is the following identity, proved in [10]; 
for completeness, the proof is also given in the Appendix. Let X be any Killing field on (Z, 7) 
and let r be any smooth symmetric bilinear form on X, satisfying the constraint equations (2.26). 
(More precisely, only the divergence-free condition in (2.26) is needed). If /i( ) is any variation of 7 
of compact support on I, then 

(3.1) J^Cxr, h {0) )dV = -2 J {S't, X)dV, 

where as in (2.29), 6' is the variation of the divergence 6 = <5 7 in the direction /i( ). (A simple 
modification of (3.1) also holds for conformal Killing fields, cf. (A. 6)). 

The relation (3.1) holds in particular for r = T( n ), where T( n ) is the stress-energy tensor associated 
to a solution g € S. Further, by Proposition 2.2 or Lemma 2.3, the equation 

(3.2) 6\ n) = -6a, 

is always solvable, for some symmetric bilinear form a on I. Clearly a is determined only up to a 
divergence-free symmetric bilinear form. 

Choose X = Z a timelike Killing field on (I, 7) and suppose the variation /i( ) of the boundary 
metric 7 has compact support; let 9S ± be any two cuts of X which enclose supp /i( )- Suppose 
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also that a = <7( n ) arises from a (global) solution h to the linearized Einstein equations, i.e. <j( n ) = 
r (n) = \n) ~ r ( ra )i ( c f- (2-29)). Then (3.1) and the divergence theorem applied to (3.2) give 

(3.3) [ {C Z T {n) ,h (0) )dV = -2 / a {n) {Z,v)dV = -2[dm+(h) - drnT (h)}, 

JX J9S± 

where dm^ is the variation of the holographic mass (2.27) at dTi^, cf. also (2.31). It is clear that 
this formula holds for any formal solution hf = (/i( ) , cr) of the linearized constraint equations (2.32) 
or (3.2): 

(3.4) f (£zr (n) , h (0) )dV = -2 f a(Z, v)dV = -2[dm + (h f ) - dm~{h f )\. 

JX J9S± 

Proof of Theorem 1.1. 

Related to the various definitions of conserved quantities for asymptotically AdS and asymptot- 
ically locally AdS spacetimes, there has been much recent discussion in the literature concerning 
the validity of the first law of black hole mechanics for such spacetimes, see in particular [34] and 
references therein and thereto. The ambiguities regarding the first law for the holographic mass 
(2.27) obtained by holographic or covariant renormalization have recently been resolved in a very 
clear analysis by Papadimitriou-Skenderis [19]. In particular, in the context of the present work, it 
is proved in [19] that if gs € £ is globally stationary, then g is a critical point for the holographic 
mass (2.27), under infinitesimal Einstein variations h with fixed boundary metric, i.e. 

(3.5) 5m = 0, 
among infinitesimal Einstein variations h such that 

(3.6) h {0) = 0. 

This is a special case of the first law of black hole mechanics in the AdS setting, namely in the case 
where there is no black hole, in the sense that the bifurcate Killing horizon is empty. 

We first observe that (3.5) also holds for all formal infinitesimal Einstein deformations hf, i.e. 
formal series solutions of the linearized Einstein equations determined by (h^,a), for a satisfying 
(2.32) on (M,g s ). To see this, from (3.3) one has: 

(3.7) J (C Z T (n) ,h (0) )dV = -2 J (6a (n) ,X)dV = -2[5m dJ:+ (h) - 5m ds -(h)}. 

Here h is any global linearized Einstein deformation with supp ht \ contained in the region of 1 
between <9£~ and <9E + and cr( n ) = fy n ) — r'^ n y Since fy ) = on <9£~, by (3.5) one has 

Sm dJ: -(h) = 0. 

Also, since (M,gs) is stationary, the left side of (3.7) vanishes. On the other hand, (3.7) holds with 
°~(n) replaced by any a satisfying (2.32), so that 

(3.8) 8mg^(h f ) = 0, 

for any formal variation hf determined by (h^,a). 

Now suppose the solution g G £ is asymptotically stationary, i.e. as U — > oo, (or t; L — > — oo), gt t 
converges to a stationary solution g^. By the discussion above, on the limit gs, one has 



Sms{hf) = / o~(Z,v)=0 
J BY. 



<9£ 
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for any formal solution hf = (h^,a) € J 7 to the linearized Einstein equations at I with fy ) = 
on <9£, (or of compact support). It follows that on the original spacetime (M,g), for \T\ sufficiently 
large, one has 

(3.9) \Sm T (h)\ = | / a(Z,u)\ < e, 

for all variations (h^,a) which are uniformly bounded by a fixed constant, independent of T and 
for which supp /i( ) is a fixed compact set in X. By Lemma 2.3, there exist such uniformly bounded 
a, for all such variations /i( ). Applying (3.4) once more to (M, g) with the data hf = (h^,a) gives 

(3.10) j^C Z T (n) , h (0) )dV = -2 jf <<5<r, Z)dV = -2[5m d ^ T+ (h f ) - 5m dl ; T _ (Mi- 
Letting T + — > +oo and T _ — > — oo, using the boundedness condition, it follows from (3.9)-(3.10) 
that 

(3.11) J(C Z T {n) ,h (0) )dV = 0, 

for all /i( ) of compact support. Since /i( ) is an arbitrary variation of the boundary metric 7 in 
I[o,t], it follows that 

(3.12) £z9(o) = and Czt^ = on X. 

Thus, the Cauchy data (g( ) , T( n ) ) or (<7( ) , <?( n ) ) for g onX are invariant under the flow of the Killing 
field Z on X. 

We claim that (3.12) together with the unique continuation property at X implies that (M,g) is 
stationary. Intuitively, this is quite clear, but the proof requires some details. 

To begin, extend Z to a vector field in the bulk by the following two-step process. First, extend 
Z to a neighborhood W of X in M by requiring 

[Zi,Vp] =0, 

with Zi|j = Z. Thus Z is extended by the flow of Vp to a vector field Zi, defined in the region 
where Vp is smooth. The corresponding form S* Z\ is then an infinitesimal Einstein deformation 
which preserves the defining function p. If <j) s denotes the flow of Z\ and g s = <p*g denotes the corre- 
sponding curve of Einstein metrics, then the geodesic defining function for g s with fixed boundary 
metric 7 is the fixed function p. Each metric g s thus has the expansion (2.9)-(2.10) with fixed 
p. Since 7 is fixed, apriori only the gi n ),s terms can vary, but (3.12) shows the variation of g^ n ),s 
vanishes to 1st order in s at s = 0. Hence, since the formal series (2.9)-(2.10) are determined by 
the (<7(o) > <7(n) ) terms, it follows that 

(3.13) (TZi = O(pP), for all p < 00, 

formally near X. More precisely, (3.13) holds to the extent that the metric g in (2.3) is smooth, (in 
a polyhomogeneous sense), up to X. In any case, one has 

(3.14) 5*Z 1 = o(p n ). 

Extend Z\ outside W arbitrarily but smoothly to all of M. 

Next we need to bring 8*Z\ into the transverse-gauge (2.22). To do this, (3.14) implies that 

uj = 55* Zi = o{p n ). 

Set Z = Z\—Y ', where Y is chosen to be the unique solution of the initial boundary value problem 



(3.15) 



55*Y = uj, 
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with zero Cauchy and boundary data; Y\s = Vo t Y\-£ = and Y\% = 0. As discussed in §2, 
this equation has a unique smooth solution at least within the domain of dependence D (£). The 
smallest indicial root of the operator 55* = ^D*D — Ric is at least (n + 1), i.e. the formal expansion 
of a solution of (3.15) has determined coefficients up to order [n + 1), cf. [35]. Since u> = o(p n ), it 
follows that 

5*Y = o(p n ). 

Hence the vector field Z on £>(£) satisfies 

(3.16) 55* Z = 0, 
with Z\q-£ the given Killing field Z on X and 

(3.17) 5*Z = o{p n ). 

Since 5*Z is an infinitesimal Einstein deformation, it follows from (2.17) that 5* Z satisfies the 
equations (2.22), i.e. 

(3.18) D*D(5*Z) - 2R(5*Z) = 0. 

Since (3.17) holds, the unique continuation property at X for (3.18) implies that 

(3.19) 5*Z = 0, 

in a neighborhood U of dT, in D (£). 

To show that Z extends to a Killing field on all of £>(£), let U s = {x G D(E) : p(x) > s > 0} 
and let C s = dU s , so that C s is a timelike cylinder. For s sufficiently small, one has 5* Z = to 
infinite order on C s and 5* Z satisfies (3.18) throughout D (£). The equations (3.18) are a hyperbolic 
system of PDEs, and at leading order are a diagonal system of scalar wave equations for which the 
boundary C s is strictly pseudoconvex. A unique continuation result of Tataru [36], then implies 
that 5* Z = in a neighborhood of C s within U s . One may then iterate this process a finite number 
of times to cover a neighborhood of the initial surface S. Of course one must change the distance 
function p near regions where p becomes singular and use instead smooth distance functions, but 
the arguments are otherwise the same. 

Since 5* Z thus vanishes to 1st order on all of £ and vanishes on X, it follows from uniqueness of 
the initial boundary value problem for (3.18) that 5* Z = on all of M. 

Observe that Z cannot become null anywhere in M. For if Z is null at some point p £ M, then 
the Killing equation (3.19) implies that flow line a of Z through p is a null line, i.e. Z remains 
null along a. Since such null lines must intersect X, this implies Z is null somewhere on X, which 
is impossible. Thus Z is timelike throughout M, so that (M, g) is globally stationary. 

■ 

Proof of Corollary 1.2. 

Theorem 1.1 applies to any time-like Killing field on (X, 7). Let /C be the cone of time- like Killing 
fields on (X, 7). Taking linear combinations, the cone K. generates the full space of Killing fields 
on X. Consider the space Sjq of solutions in £ which are invariant under an effective /C-action, 
restricting to the action of /C at conformal infinity. Then Theorem 1.1 implies that any solution 
(M, g) G £ which to the future and past is weakly asymptotic to (possibly distinct) elements in Sx, 
is necessarily a fixed solution in Sic- 

Applying this to metrics in £ whose conformal infinity is the Einstein static cylinder, it follows 
that the component of the identity of the isometry group of the Einstein static cylinder IR x S"" -1 
extends to a group of isometries of any (M,g) G £ . In particular, any such (M,g) has an isometric 
R x SO(n) action. This implies that the Einstein equations (1.1) reduce to a system of ODE's and 
it is standard that the only globally smooth solution of this system is the exact AdS solution. 

■ 
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Exactly the same arguments can be applied to spacetimes (M, g) £ £ whose conformal infinity 
is homogeneous, i.e. (X, 7) has a transitive group of isometries. The Einstein equations for (M,g) 
then reduce to a system of ODE's, (in the variable p), and the requirement that the solutions 
are smooth in the interior typically gives either a unique solution, or uniqueness up to a set of 
parameters which determine the topology of (M,g). 

To illustrate on some concrete examples, consider the AdS soliton metric of Horowitz- Myers [11]. 
In the toroidal compactification, (X, 7) is the flat product metric on R x T n_1 on the (n — l)-torus. 
This is clearly homogeneous. The corresponding ODE's may be solved explicitly and have a unique 
smooth solution on (M,g) ~ R x D 2 x T n ~ l £ £ up to the choice of topology, (a choice of the 
disc D 2 bounding an S 1 C T n_1 ). This proves uniqueness of the AdS soliton metric among all 
(dynamical) metrics in £ with the given conformal infinity and topology which are asymptotic at 
t = +00 and t = —00 to an AdS soliton metric. (This argument can be extended without difficulty 
to the case where X is compactified to a single circle instead of the full (n — l)-torus). Exactly the 
same results hold for the recent AdS soliton metric analysed by Copsey- Horowitz [12]. 

Similar uniqueness results also hold with respect to perturbations of such homogeneous conformal 
infinities. For example, suppose (M, g) is an asymptotically simple, globally static solution of the 
Einstein equations (1.1) with conformal infinity (X, 7) which is non-degenerate, (e.g. (M,g) has 
non-positive curvature). Given a static or stationary perturbation (X, 7) of the boundary data 
(X, 7), there is a unique globally static or stationary asymptotically simple solution (M,g) close to 
(M,g) with conformal infinity (X, 7), cf. [37], [38]. Theorem 1.1 then implies the solution (M, g) is 
unique among all dynamical solutions to the Einstein equations in £ which are asymptotic to the 
future and past to the given static or stationary solution. 

Remark 3.1. It is worth emphasizing that the results above require the solutions (M,g) of the 
Einstein equations to be globally defined and non-singular. On the one hand, this is apparent 
from the proof. Theorem 1.1 uses the global vanishing (3.5) of the variation of mass on stationary 
spacetimes; per se, this is false if there are inner boundaries in addition to the boundary at conformal 
infinity. 

In fact, Theorem 1.1 or Corollary 1.2 are false for solutions of the Einstein equations defined 
only in a neighborhood of conformal infinity (X, 7). For example, let (X, 70) be the Einstein static 
cylinder and let T( n ) be any analytic symmetric bilinear form on X which is asymptotic to as 
t —> ±00, and which satisfies the constraint equations (2.26). The result of Kichenassamy [22] 
mentioned above gives the existence of a solution of the Einstein equations (1.1) defined in a 
neighborhood of X with the given T( n ), i.e. <7( n ) term, on X. This solution is asymptotic to the exact 
AdS solution at t = ±00, but is not exact AdS unless T( n ) = 0. 



4. Discussion 

In the context of the Euclidean (or Riemannian) version of the AdS/CFT correspondence, it is 
important to know to what extent the boundary data (dM, 7) determine the bulk solutions (M, g) 
of the Einstein equations (1.1). Although it is possible in general that there are infinitely many 
topological types for (M,g), or that, fixing the topology, the space of solutions has infinitely many 
components, there is only at most a finite dimensional moduli space of solutions when one fixes 
the topology and component. This follows essentially from the elliptic character of the Einstein 
equations (1.1). Thus, the "Cauchy data" (<?(o) , <7( n ) ) uniquely determine the bulk solution (M,g) 
up to local isometry, (cf. [10]), and given gt \, although the stress-energy term gr n \ may not quite 
be uniquely determined by the boundary metric g( ), it is determined up to a finite dimensional 
space of parameters, (given a choice of topology and deformation component). 
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In this Euclidean context, the fact that isometries of the boundary (dM, 7) necessarily extend to 
isometries of any smooth global bulk solution (M, g) is a simple and clear illustration of (elementary 
or classical) aspects of the AdS/CFT correspondence. 

On the other hand, in the context of the Lorentzian solutions of the Einstein equations (1.1) in 
£ , there is an infinite dimensional space of normalizable modes, i.e. L 2 solutions of the linearized 
Einstein equations. At least in some situations, for example that of exact AdS spacetime, these 
linearized solutions remain uniformly bounded in time, i.e. gAdS 1S linearization stable. One expects 
that that gAdS 1S i n f &c t dynamically stable, with the behavior of the nonlinear exact solutions 
nearby to gAdS well-modeled on the linearized behavior. 

In fact, Friedrich [39] has shown that one may solve the initial boundary value problem for the 
Einstein equations (1.1) at least locally in time, in 3+1 dimensions. Thus, suppose the boundary 
data (1, 7) are globally stationary, and that £ is a Cauchy surface with initial data (go , K) satisfying 
the constraint equations (2.15) and matching 7 at the corner <9£ C I, (cf. also [40]). Then there is 
a solution g of the Einstein equations (1.1), with conformal infinity (2, 7) defined on a thickening 
J x E of E, realizing the given Cauchy data. This gives an infinite dimensional space of local-in- 
time exact solutions. If the data (go,K) are (arbitrarily) close to the data (g_i,0) of exact AdS 
spacetime, where g-\ is the hyperbolic metric on the n-ball, then the resulting solution will exist for 
an (arbitrarily) long time, to the future and past. One would expect, although this remains to be 
proved, that such solutions extend to global-in-time solutions which remain globally asymptotically 
simple and globally close to gAdS- 

Theorem 1.1 implies that such global solutions cannot be both future and past asymptotic to 
the exact AdS spacetime. This is of course in strong contrast to the asymptotically flat situation 
(A = 0), where the Christodoulou-Klainerman theorem [2] implies that small global perturbations 
of Minkowski spacetime disperse in time and tend to the flat solution, preserving, (at least to a 
certain degree), asymptotic simplicity and boundedness; in some situations asymptotic simplicity 
and boundedness is preserved to all orders, [5]. To a certain extent, it is the global conformal 
or causal structure which leads to these differences. Thus, observe that the Bondi mass on Z + is 
not preserved in time, but decreases monotonically. While the ADM mass is preserved under time 
evolution, it is defined with respect to the singular structure at spacelike infinity t°, which does not 
match smoothly with the geometry of null infinity X + . 

Similarly, with regard to the global stability results of de Sitter spacetimes in [8] and [9] , here it 
is not clear if there is even a reasonable definition of mass, since X is spacelike. Given a definition of 
mass, as for instance in [41], it is not directly related to the dynamical evolution of the spacetime; 
it is a fixed quantity at future or past spatial infinity. 

On the other hand, Theorem 1.1 is consistent with the known oscillatory behavior of the Einstein 
equations linearized at exact AdS, cf. [1]. It is unknown if general asympotically locally AdS Einstein 
spacetimes are linearization stable, or what conditions will guarantee linearization stability. 

To conclude, we discuss informally the implication of Theorem 1.1 on the global dynamical 
behavior of bounded solutions in 8. Thus, consider the subspace £k of uniformly bounded solutions 
in £, i.e. solutions which remain uniformly bounded by a fixed constant K in time, in a suitable 
smooth norm. Then Sk is weakly compact, in that in any sequence in £k has a subsequence 
converging in a weaker topology to a limit solution again in £x- Standard results in dynamical 
systems, (cf. [42] for example), then imply that the resulting Einstein flow on the constraint space 
Ck has (many) uniformly recurrent points: thus there exist times U — > 00, (or ti — > —00), with 
U+i — ti<T, for some fixed T, and points (g^ , K ti ) G Ck such that 

dist((g u+1 ,K u+1 ),(g u ,K u )) < e, 
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where the distance is taken in the weaker norm. Here e may be chosen arbitrarily small, but then 
with T becoming possibly arbitrarily large. 

Theorem 1.1 implies that no orbits of the flow on Ck are asymptotic to fixed points as t — > ±00. 
This indicates that most all orbits in Ck should be almost periodic. 

It is an open question whether there exist periodic (non-stationary) solutions in £. The proof 
of Theorem 1.1 can be easily adapted to rule out periodic solutions provided, for any periodic 
boundary variation h^, there exists a periodic solution a of the linearized constraints as following 
(3.7). However, it is not clear that such periodic a always exist; I am very grateful to Bob Wald 
for discussions on this point. 

Appendix 

Let I be a Killing field with respect to a Lorentz metric (X, 7), and let r be a divergence- free 
symmetric bilinear form on T, 5t = 0. We give here a proof of the identity (3.1) from [10], i.e. 

(A.l) f (C X T, K )dV = -2 f (S'(r),X)dV, 

Jx Jx 

where k = ■^('j+uk)\ s= o is a variation of 7 of compact support and 5' = -^5 1+sk \ s= q is the variation 
of the divergence. The proof of (A.l) below holds for metrics of any signature. 
To prove (A.l), we use the following standard formulas, cf. [43] for example: 

(A.2) C v <f) = V V (j) + 2W o 0, 



(A.3) {5*)'V = \V v k + 5*V ok, 

for any vector field V. Here <fi o tp is the symmetrized product; (4> o tp)^ = |((<fo, tpj) + (<fij, ipi)) and 
(6*V) ij = ±(V itj + V jti ). 

To begin, by (A.2), f x (£ X T,K) = /^Vjt, k) + 2(VIot,k). Since k is symmetric, (VXor, n) = 
(5*Xot, k). For the first term, write (Vxt, k) = X(t, k) — (r, Vxt). The first term here integrates 
to SX(t, k), (here we use the fact that k is of compact support), while by (A.3), the second term 
is -(r, V x n) = -2(r, (S*)'X) + 2(r, 5*X o k). Hence 

f (C x t, n)dV = -2 f (T,(5*)'X)dV + 4 f (t,5*Xok)+ j 5X{r,K)dV. 

Next, a straightforward computation using the fact that 8t = gives 

J (r, (5*)'X)dV = J ((<5')(r), X)dV + 2 J [{to 5*X, k) - i(r, 5*X)trK]dV : 

the last two terms come from variation of the metric and volume form. Combining these computa- 
tions gives 

(A.4) j {£ x t, n)dV = -2 j (S'(T),X)dV + j [SX{t, k) + (r, 5*X)trK]dV. 

This gives (A.l) when X is Killing, i.e. 5*X = 0. 

Let txl denote the conformal Killing operator: £x7 = £x7 — 2d ™ x " 7, (where divX = —SX). 
Then the calculations above give, (again for divergence- free r), 

(A.5) f {C x t + divXr, n)dV = -2 f (S'(t), X)dV + \ I (t, Cxl)dV + \ j divX a ■ tr ndV, 
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where a = trr. If a satisfies the linearized constraint equations (2.32) then a simple calculation 
from (A. 5), gives 

(A.6) J (C x t + [1- l]divXT,n)dV = 

j (a + \trKT 1 Cxl)dV -2 J div(a{X))dV + ± J divX(a ■ tr k + 2a')dV. 

When X is a conformal Killing field on (X, 7), the first term on the right vanishes. When r = T( n ), 
the second term gives the difference of the variation of the mass when X is timelike by the divergence 
theorem while the integrand on the left is related to the transformation properties of the stress- 
energy tensor under conformal changes, cf. [17]. 
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